Abstract of DE 1 0038460 (A1 ) 

Translate this text 

A fluid flow sensor uses Particle 
Image Velocimetry to plot the 
vector flow (20) of fluids such as 
the charge in an engine cylinder 
with data processing to determine 
the real center of turbulent 
rotation (38). 
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Die folgenden Angaben sind den vom Anmelder eingereichten Unterlagen entnommen 

@ Verfahren zum Messen von Stromungsverhaltnissen von Fluiden 

(57) Die Erfindung betrifft ein Verfahren zum Messen von 
Stromungsverhaltnissen von Fluiden, insbesondere einer 
Ladungsmasse in Zylindern einer Brennkraftmaschine r 
wobei die Stromung eine Wirbelstromung, insbesondere 
eine Tumb j e- und/oder eine Drallstromung ausbildet, und 
aus aufeinanderfolgenden Momentanaufnahmen der 
Stromungsverhaltnisse jeweilige Stromungsvektoren mit 
Stromungsrichtung und Stromungsgeschwindigkeit an 
verschiedenen Orten bestimmt und daraus Kennzahlen K 
der Stromung berechnet werden. Hierbei werden aus den 34 > 
Momentanaufnahmen jeweilige tatsachliche Mittelpunk- 
te der Wirbelstromung ermittelt und die Kennzahlen aus 
den Stromungsvektoren an den verschiedenen Orten be- 
rechnet, wobei zur Berechnung fur jeden Ort eines jewei- 
ligen Strornungsvektors ein zugehoriger Ortsvektor ver- 
wendet wi 'd, der vom tatsachlichen Mittelpunkt der Wir- 
■ belstrbmung als Ursprung ausgeht. 
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Beschreibung 

[0001] Die Hrfindung betrifft ein Verfahren zum Messen von Slromungsverhaltnissen von Fluiden, insbesondere einer 
Ladungsmasse in Zylindern einer Brennkraftmaschine, wobei die Stromung eine Wirbelstromung, insbesondere eine 
5 Tumble- und/oder eine Drallstromung ausbildet, und aus aufeinanderfolgenden Momentanaufnahmen der Stromungs- 
vcrhaltnissc jcweiligc Stromungsvektoren mit Stromungsrichtung und Stromungsgeschwindigkcit an vcrschicdcncn Or- 
len beslimmt und daraus Kennzahlen K der Stromung berechnet werden, gemaB dem Oberbegriff des Anspruchs 1. 
[0002] In verschiedenen technischen Gebieten ist es erstrebenswert, eine Fluidbewegung in einem abgeschlossenen 
Raum zu beobachten und zu analysieren. So ist beispielsweise die Ladungsbewegung in einem Brennraum bzw. Zylinder 
10 sowohl fur Otto- als auch Dicselmotoren ein Parameter, der die Verbrennung und damit die Wirtschaftlichkeit und Lei- 
stungsfahigkeit wesenllich beeinfluBt. Grundsatzlich wird zwischen Drall, einer Ladungsbewegung urn die Zylinder- 
achse, und Tumble, einer Wirbelstromung urn eine Achse senkrecht zur Zylinderachse, unterschieden, Meist treten beide 
Ladungsbewegungen gleichzeitig in verschiedener Form auf. 

[0003J Aus "Motortcchnische Zeitschrift" MTZ 61 (2000) 1, S, 40 ff ist cs bekannt, Drall- und Tumblckennzahlcn mit 

15 Ililfe von mechanise h arbeitenden Aufnahme verfahren zu ermitteln. Zum Erfassen beider Arten der Ladungsbewegung 
werden an einem Prufstand fur stationare Stromung integrative Methoden angewendet. Nachteilig ist jedoch, daB mecha- 
nische Hilfsmittel, wie ein Fliigelrad oder MomentenmeBgerat, zur Aufnahme von Stromungsrichtung und Stromungs- 
geschwindigkeit verwendet werden und diese in der zu messenden Stromung selbst angeordnet werden mussen, wodurch 
die Stromtng in unerwunschter Weise beeinfluBt wird. 

20 [0004] Aus der EP 0 566 120 Bl ist es bekannt, zur Messung der Stromungsverhaltnisse u. a. in Brennraumen von 
Brennkraftmaschinen das sogen. PIV- Verfahren (Particle Image Velocimetry) anzuwenden. Das stromende Fluid enthalt 
hicrbei Markierungspartikel (tracer partikels) und das Fluid wird mittels eines kontinuierlichen oder gepulsten Lasers an- 
gestrahlt. Dies machl die Fluidverteilung, Stromungsrichtung und -geschwindigkeit mittels der Markierungspartikel 
sichtbar. Entsprechende Bilder werden aufgenommen und mittels einer Bildverarbeitung ausgewertet. Es werden bei- 

25 spielsweise mit 25 Hz zwei dimension ale Bilder in Form von Vektorfeldern aufgenommen und anschlieBend anhand ei- 
ner Abfolge von mehrcrcn aufeinanderfolgenden Vektorfeldern die tatsachlichen Stromungsgeschwindigkciten ermittelt. 
Bei dem PIV- Verfahren in Verbindung mit Tumble-Messungen wird das dreidimensionale Stromungsfeld im Zylinder 
auf die zweidimensionale Ebene reduziert, welche im groBten Zylinderdurchmesser vorliegt. Fur die Ermittlung der 
Drall zahl wird ein urn 90° gedreht.es MeB verfahren angewandt, wobei das Stromungsfeld senkrecht zur Zylindermitten- 

30 achse ausgewertet wird. Hicrbei miBt beispielsweise eine Kamera senkrecht von unten in den Brennraum. 

[0005] Bei der anschlieBenden mathematischen Berechnung eines Tumble- bzw. Drallkennwertes wird jedoch bei alien 
vorgenanmen Verfahren auf einem Ortsfesten Mittelpunkt der Stromung aufsetzt, welcher aber nur theoretisch gegeben 
und im tatsachlichen Bctricb nicht an der angenommcnen Stclle vorhanden ist. Daher wcichen die mit den bckannten 
Methoden ermittelten Werte zum Teil erheblich von den real vorhandenen Werten ab. 

35 [0006] Der vorliegenden Erfindung liegt die Aufgabe zugrunde, ein verbessertes Verfahren der obengenannten Art zur 
Verftigung zu stellen, wobei die obengenannten Nachteile Liberwunden werden und eine noch genauere Messung und 
Charakterisicrung von Wirbelstromungcn crziclt wird. 

[0007J Diese Aufgabe wird erfindungsgemaB durch ein Verfahren der o. g. Art mit den in Anspruch 1 gekennzeichne- 
ten Merkmalen gelost. Vorteilhafte Ausgestaltungen der Erfindung sind in den abhangigen Anspruchen angegeben. 
40 [0008] Dazu ist es erfindungsgemaB vorgesehen, daB aus den Momentanaufnahmen jeweilige tatsachliche Mittel- 
punkte der Wirbelstromung ermittelt und die Kennzahlen aus den Stromungsvektoren an den verschiedenen Orten be- 
rechnet werden, wobei zur Berechnung fur jeden Ort eines jeweiligen Stromung svektors ein zugehbrigcr Ortsvektor ver- 
wendet wird, der vom tatsachlichen Mittelpunkt der Wirbelstromung als Ursprung ausgeht. 

[0009] Dies hat den Vorteil, daB eine Bewegung des gesamten Stromungsfeldes in der Betrachtungsebene unkritisch 
45 und ohne EinfluB auf das MeBergebnis ist. Das erfindungsgemaBe Verfahren benotigt keinen festen Mittelpunkt. 

[0010] Urn die zu messende Stromung selbst so wenig wie moglich zu beeinflussen ist cs besonders vorteilhaft, die 
Momentanaufnahmen mittels eines PIY-Verfahrens (Particle Image Velocimetry) herzustellen. 
[0011] Ei ne Kennzahl der Wirbelstromung wird gemaB folgender Formel berechnet: 
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L: Langc cles Zylindcrs, 

D: Dure hn lesser des Zylinders, 

r\\ Abstand zwischen MeBpunkt i und realer Achse der Wirbelstromung, 

l x \ Ortsvektor des MeBpunktes i ausgehend von der realen Achse O der Wirbelstromung, 

(O exp : Winkclgcschwindigkcit der gesamten McBflache, 30 
<j»moc* Winkelgeschwindigkeit des Motors, 
Q: DurchfluBwert, 

V m ean: mitilere Geschwindigkeit, folgend aus dem DurchfluBwert Q und Zylinderdurchmesser D, 
\{. Geschwindigkeitsvektor am MeBpunkt i, 

n: Anzahl von MeBpunkten in der MeBflache (Vektorfeld), 35 
N: Anzahl von gemessenen MeBflachen (Vektorfeldern), 
K: Kennzahl der Wirbelstromung. 

[0012] Als Kennzahl K wird eine Tumblezahl T, welche eine Wirbelstromung in einem Zylinder einer Brennkraftma- 
schine mit Drehachse senkrecht zur Zylinderachse beschreibt, und/oder eine Drallzahl DZ, welche eine Wirbelstromung 
in einem Zylinder einer Brennkraftmaschine mit Drehachse parallel zur Zylinderachse beschreibt, berechnet. 40 
[0013] Durch das crfindungsgcmaBc Vcrfahrcn ist die Berechnung einer zusatzlichcn charaktcristischen GroBe bzw, 
Kennzahl moglich, namlich die Bcstimmung einer Stabilitatszahl S, welche die Bewcgung des gesamten Stromungsfcl- 
des beschreibt und gemaB folgender Formel berechnet wird: 

S = 45 
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und 

r^: Abstanc zwischen realer Wirbelachse i und idealer Wirbelachse O, 65 

rroeai)' milUercr Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 

r^v! Abweichung Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 

r raax : maximale Abweichung des Abstands zwischen realer Wirbelachse i und idealer Wirbelachse O, 
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D: Durchmesser des Zylinders 

N: Anzahl von gemessenen MeBflachcn (Vektorfcldern) 

S: Stabilitalszahl. 

[0014] Weitere Merkmale, Vorteile und vorteilhafte Ausgestaltungen der Erfindung eigeben sich aus den abhangigen 
5 Anspruche i, sowie aus der nachstehenden Beschreibung der Erfindung anhand der beigefugten Zeichnungen, Diese zei- 
gen in 

[0015] Fig. 1 eine schematische Darstellung eines stationaren Prufstandes fiir einen Zylinderkopf, 

[0016] Fig. 2 eine graphische Veranschaulichung des MeBprinzips des PIV-Verf aniens (Particle Image Velocimetry), 

[0017] Fig. 3 eine beispielhafte Momentaufnahme eines MeBfeld zur Berechnung einer Tumblezahl, 

10 [0018] Fig, 4 cine beispielhafte Momentaufnahme eines MeBfeld zur Berechnung einer Drallzahl und 
[0019] Fig. 5 eine beispielhafte Momentaufnahme eines MeBfeld zur Berechnung einer Stabilitatszahl. 
[0020] Nachfolgend wird die Erfindung anhand der Ermittlung von Kennzahlen einer Ladungsmasse in einem Kolben 
einer Brennkraftmaschine naher beschrieben. Dies ist jedoch lediglich beispielhaft zu verstehen. Das erfindungsgemaBe 
Verfahren ist zur Analyse jeglicher Stromungen gecignet, die eine Wirbelstromung ausbilden, 

15 [0021] Die Fig. 1 und 2 veranschaulichen das PIV- Verfahren (Particle Image Velocimetry). Wie insbesondere aus Fig. 
1 ersichtlich, wird ein Zylinderkopf 10 mit EinlaBkanalen 12 auf einem Prufstand montiert, der einen Zylinder 14 mit ei- 
nem feststehenden flachbodigen Kolben (nicht dargestellt) ausbildet. Bei einem konstanten Druckunterschied zwischen 
dem Eintriit der EinlaBkanalc 12 und dem Zylinder 14 wird cine Kanalstromung in den Zylinder 14 eingeleitet, Mittels 
einer Kamera 16 wird mil einem PIV- Verfahren, welches nachfolgend unter Bezugnahme auf Fig. 2 naher erlautert wird, 

20 ein zweidimensionales MeBfeld 20 aufgenommen. In einem Computer 18 werden an verschiedenen Orten dieses MeB- 
feldes die entsprechenden Geschwindigkeitsvektoren, d. h. Betrag und Richtung der Geschwindigkeit, analysiert. 
[0022] Wie aus Fig. 2 ersichtlich, wird mittels des PIV-Vcrfahrcns das sich unterhalb des Zylindcrkopfes 10 ausbil- 
dende wSlrdmungsfeld 26 gemessen. Ilierzu strahlt ein Laser 22 durch eine Lichlschniltoptik 24 auf das Stromungsfeld 26, 
so daB sich ein Lichtschnitt 28 ergibt, der das Stromungsfeld in einer zweidimensionalen Ebene 30, d. h. in einem Beob- 

25 achtungsfeld 30, fiir die Kamera 16 sichtbar macht. Die Kamera 16 nimmt dabei das Beobachtungsfeld 30 durch ein Ob- 
jektiv 32 auf. Fiir mehrere Zeiten t wird auf diese Weise ein zweidimensionales MeBfeld 20 aufgenommen. Entspre- 
chende Tumble- und Drallzahlen werden in nachfolgend noch naher beschriebener Weise durch numerische Integration 
der fiir verschiedene Zeiten t gemessenen zweidimensionale MeBfeld 20 aus den Geschwindigkeitsvektoren und deren 
Lage berechnet. 

30 [0023] Fig. 3 vcranschaulicht die Berechnung einer TYimblczahl, welche eine Wirbelstromung urn cine Achse senk- 
recht zu einer Langsachse des Zylinders 14 beschreibi. Der Zylinder weist eine Lange L 34 sowie einen Durchmesser D 
36 auf. Hierzu wird der Lichtschnitt 28 gemaB Fig. 2 in einer Ebene parallel zur Langsachse des Zylinders 14 ausgerich- 
tct und zwar am On des groBtcn Durchmessers des Zylinders 14, d. h. die Zylinderlangsachse verlauft genau in der 
Ebene des Lichtschnittes 28. Hierdurch wird die Tumblestromung im Zylinder 12 erfaBt, die eine Wirbelstromung urn 

35 eine reale Achse O 38 darstellt, die eine Winkelgeschwindigkeit co exp (t) aufweist. Eur mehrere Zeiten t werden entspre- 
chende zweidimensionale MeBfelder 20 aufgenommen. An verschiedenen Orten i mit Ortsvektor r j(t.) undendlicherFla- 
che A, wird ein entsprechender GeschWindigkeitsvektor V^t) bestimmt. Aus diesen Geschwindigkeitsvektoren Vi(t) wird 
dann eine Winkelgeschwindigkeit COj am MeBpunkt i ermittelt. Aus alien Winkelgeschwindigkeiten (Oi an den verschie- 
denen Mefcpunkten i ergibt sich dann die Winkelgeschwindigkeit a> exp (t) der gesamten MeBflache. Das Verhaltnis von 

40 co cxp (t) zu einer Winkelgeschwindigkeit (0 raot der Brennkraftmaschine ergibt dann die Tumblezahl. 
[0024] Die Berechnung der Tumblezahl erfolgt nach folgenden Formcln: 



v. 



Q 



mean 



ttD 2 /4 



45 



V 



CO 

i 



mean 



mot 



2L 



50 



n 



55 




60 




wobei 



65 L: Lange des Zylinders, 

D: Durchmesser des Zylinders, 

T\\ Abstand zwischen MeBpunkt i und realer Achse der Wirbelstromung, 

7 4 : Qrtsvektor des MeBpunktes i ausgehend von der realen Achse O der Wirbelstromung, 
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G) exp : Winkelgeschwindigkeit der gesamten MeBflache, 
CO raot : Winkelgeschwindigkeit des Motors, 
Q: DurchfluBwert, 

Vmean^ mittlere Geschwindigkeit, folgend aus dem DurchfluBwert Q und Zylinderdurchmesser D, 

V;: Geschwindigkeitsvekror am MeBpunkt i, 5 
n: Anzahl von McBpunkten in der MeBflache (Vektorfcld), 
N: Anzahl von gemessenen MeBflachen (VekLorfeldern), 
T: Tumble. 

[0025] Das auBergewohnliche an dem erfindungsgemaBe Verfahren liegt darin, daB im Gegensatz zu herkommlichen 
Bcrechnungsvcrfahrcn nicht cine vorbestimmte, feststchende idealc Drehachse als Ursprung fur die Ortsvektoren rtft) 10 
verwendel wird, sondern slatl dessen derjenige Punkt O 38 in dem jeweiligen zweidimensionalen MeBfeld 20, bei dem 
die reale Drehachse der Tumblestromung das zweidimensionale MeBfeld 20 schneidet. 

[0026] In analoger Weise wird die Drallzahl DZ gemaB der Darstellung in Fig. 4 ermittelt, so daB zur naheren Rrlaute- 
rung auf die obige Beschrcibung der Fig. 3 verwiesen wird. Im Unterschied zur Krmittlung der 'I\imblezahl T gemaB Fig. 
3 ist der Lxhtschnitt senkrecht zur Langsachse des Zylinders 14 ausgerichtet, so daB eine Wirbelstromung (Drall) des 15 
Fluides um eine Achse aufgenommen wird, die parallel zur Zylinderlangsachse verlauft. 
[0027] Die Berechnung der Drallzahl erfolgt nach folgenden Formeln: 
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wobei 

L: Lange des Zylinders, 

D: Durchrr esser des Zylinders, 

r{. A b stand zwischen MeBpunkt i und realer Achse der Wirbelstromung, 45 
l x \ Ortsvektor des McBpunktes i ausgehend von der rcalen Achse O der Wirbelstromung, 
«>exp' Winkelgeschwindigkeit der gesamten MeBflache, 
(o mo[ : Winkelgeschwindigkeit des Motors, 
Q: DurchfljBwert, 

V mcan : mittlere Geschwindigkeit, folgend aus dem DurchfluBwert Q und Zylinderdurchmesser D, 50 

V\\ Geschwindigkeitsvektor am MeBpunkt i, 

n: Anzahl von MeBpunkten in der MeBflache (Vektorfeld), 

N: Anzahl von gemessenen MeBflachen (Vektorfeldcrn), 

DZ: Drallzahl. 

[0028] Wiederum wird als Ursprung fur die Ortsvektoren fi(t) derjenige Punkt O 38 in dem jeweiligen zweidimensio- 55 
nalen MeBfeld 20 gewahlt, bei dem die reale Drehachse der Drallstromung das zweidimensionale MeBfeld 20 schneidet. 
[0029] Z asatzlich wird durch das erfindungsgemaBe Verfahren die Berechnung ciner zusatzlichcn KcnngroBe moglich, 
namlich der Stabililaiszahl, die die Bewegung des gesamten Stromungsfeldes beschreibt. 
[0030] Die Berechnung der Stabilitatszahl S erfolgt nach folgenden Formeln: 

60 
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15 




T[. Abstand zwischen realer Wirbelachse i und idcalcr Wirbclachse O, 
fmean- mittlerer Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 
r dev : Abweichung Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 
25 r inax : maximale Abweichung des Abstands zwischen realer Wirbelachse i und idealer Wirbelachse 0, 
D: Durchmesscr des Zylinders 

N: Anzahl von gemessenen MeBfiachen (Vektorfeldem) 
S: Stabilitatszahl. 

[0031] Hierbei entspricht der Ursprung O 40 der idealen Wirbelachse. In der Fig. 5 reprasentieren die Kreuze die ver- 
30 schiedcnen rcalen Drehachsen in den zeitlich beabstandet aufgenommenen zweidimensionalcn MeBfeldern. Es wird ein 
Mittelwert aller Abweichungen der realen Wirbelachsen von der idealen Wirbelachsen gebildet. Das Verhaltnis dieses 
Mittelwert'js zu einem maximalen Ahstand ergibt die Stabilitatszahl S. 

BEZUGSZEICHENLISTE 

35 

lOZylinderkopf 
12 EinlaBkanal 
14 Zylinder 
16 Kamera 
40 18Compuier 

20 zwcidinicnsionales McBfcld 
22 Laser 

24 Lichtschnittoptik 

26 Stromungsfeld 
45 28 Lichtschnitt 

30 zweidimensionalcn Ebcne/Bcobachtungsfcld 

32 Objektiv 

34 Lange L 

36 Durchmesscr D 
50 38 reale Achse der TLimblestiomung 

40 Ursprung fur Stabilitatszahl 

DZDrallzahl 

K Kennzahl 

TTumblezahl 

55 

Patentanspriiche 

1. Verfahren zum Messen von Stromungsverhaltnissen von Fluiden, insbesondere einer Ladungsmasse in Zylin- 
dern einer Brennkraftmaschine, wobei die Stromung eine Wirbelstromung, insbesondere eine Tumble- und/oder 

60 eine Drallstromung ausbildet, und aus aufeinanderfolgenden Momentanaufnahmen der Stromungsverhaltnisse je- 

weilige Stromungsvektoren mit Stromungsrichtung und Stromungsgcschwindigkeit an verschiedcncn Orten be- 
stimnit und daraus Kennzahlen K der Stromung berechnet werden, dadurch gekennzeichnet, da!3 aus den Momen- 
tanaufnahmen jeweilige tatsachliche Mittelpunkte der Wirbelstromung ermittelt und die Kennzahlen aus den Stro- 
mungsvektoren an den verschiedenen Orten berechnet werden, wobei zur Berechnung fur jeden Ort eines jeweili- 

65 gen Slromungsvektors ein zugehoriger Ortsvektor verwendet wird, der vom tatsachlichen Mittelpunkt der Wirbel- 

slromung als Ursprung ausgeht, 

2. Veifahren nach Anspruch 1, dadurch gekennzeichnet, daB die Momentanaufnahmen mittels eines PlV-Verfah- 
rens (Particle Image Velocimetry) hergestellt werden. 
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3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, daB eine Kennzahl der Wirbelstromung gemaB fol- 
gcnde r Formel berechnet wird, 

l-\ ( °mot 



wobei 
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und 
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und 

Vmean= 'xD 2 /4 

' 25 

und 

L: Lange des Zylinders, 

D: Durchmesser des Zylinders, 

IV Abstand zwischen MeBpunkt i und realer Achse der Wirbelstromung, 

x\\ Orsvektor des McBpunktes i ausgehend von der realcn Achse O der Wirbelstromung, 30 
(De* P : Winkelgeschwindigkeit der gesamten MeBflache, 
(0, noc : Winkelgeschwindigkeit. des Motors, 
Q: DurchfluBwcrt, 

V incan : mittlere Geschwindigkeit, folgend aus dem DurchfluBwert Q und Zylinderdurchmesser D, 

vV Geschwindigkeitsvektor am MeBpunkt i, 35 

n: Anzahl von MeBpunkten in der MeBflache (Vektorfeld), 

N: Anzahl von gcmessenen McBflachen (Vcktorfeldern), 

K: Kennzahl der Wirbelstromung. 

4. Veifahren nach einem der vorhergehenden Anspriiche, dadurch gekennzeichnet, daB als Kennzahl K eine 1\im- 
blezabl T, welche eine Wirbelstromung in einem Zylinder einer Brennkraftmaschine mit Drehachse senkrecht zur 40 
Zylincterachse beschreibt, und/oder eine Drallzahl DZ, welche eine Wirbelstromung in einem Zylinder einer Brenn- 
kraftmaschine mit Drehachse parallel zur Zylinderachse beschreibt, berechnet wird. 

5. Veifahren nach einem der vorhergehenden Anspruche, dadurch gekennzeichnet, daB zusatzlich eine Stabilitats- 
zahl S berechnet wird gemaB der Formel, 

r 45 

r max 
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50 
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und 

T['. Abstand zwischen realer Wirbelachse i und idealer Wirbelachse O, 
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r rtl eajp mittlerer Abstand /wischen realer Wirbelachse i und idealer Wirbelachse 0, 

r cev : Abwcichung Abstand zwischenrcalcr Wirbelachse i und idealer Wirbelachse 0, 

rn, ax : maxiniale Abweichung des Abstands zwischen realer Wirbelachse i und idealer Wirbelachse 0, 

D: Durchmesser des Zylinders 

N: Anzahl von gemessenen MeBflachen (Vektorfeldern) 
S: S tsibili tatszahl. 
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